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Electric field effect devices based on mesoscopic graphite are fabricated for galvanomagnetic mea- 
surements. Strong modulation of magneto-resistance and Hall resistance as a function of gate voltage 
is observed as sample thickness approaches the screening length. Electric field dependent Landau 
level formation is detected from Shubnikov de Haas oscillations in magneto-resistance. The effective 
mass of electron and hole carriers has been measured from the temperature dependant behavior of 
these oscillations. 



Graphite is a semimetal with highly anisotropic elec- 
tronic structure featuring nearly compensated low den- 
sity electrons and holes with very small effective mass Q . 
Such an unusual electronic structure is the basis of unique 
electronic properties of other graphitic materials, such as 
fullerenes and carbon nanotubes 2], and may lead to 
novel manifestations in two-dimensional graphene mate- 
rials. For this reason, electron transport in graphite has 
recently been the subject of extensive theoretical 0, 
and experimental 0, S S 0] investigations. In par- 
ticular, interesting size dependent galvanomagnetic ef- 
fect has been observed [lfj, yjj in thin layers of graphite, 
with thickness approaching ~ 10 nm. On this mesoscopic 
length scale, the electrostatic field-effect (EFE) modula- 
tion of the charge carrier concentration is expected to be 
very effective, owing to the low density of nearly com- 
pensated carriers in graphite. However, the EFE depen- 
dent galvanomagnetic measurements in graphite have not 
been carried out in previous studies due to the difficulty 
in obtaining adequate sample geometries. 

In this letter, we present results from the magne- 
toresistance (MR) and Hall resistance measurements in 
mesoscopic graphite crystallites consisting of as few as 
~ 35 atomic layers. Strong modulation of galvanomag- 
netic transport has been observed as the gate electric 
field changes. EFE dependent Shubnikov de Hass (SdH) 
oscillations, signatures of Landau level formation of elec- 
trons and holes, have been observed at low temperatures. 
In addition, the effective mass of electrons and holes are 
measured by investigating the temperature damping of 
SdH amplitudes for each type of carriers. 

The mesoscopic graphite devices used in this exper- 
iment are fabricated using a unique micro-mechanical 
method. The details of the device fabrication are de- 
scribed elsewhere (l2|. In brief, small graphite crys- 
tal blocks are extracted from bulk highly oriented py- 
rolytic graphite (HOPG) using micro-patterning followed 
by micro-mechanical manipulation. A detached HOPG 
block is then transferred and fixed onto a micro- machined 
Si cantilever. By operating an atomic force microscope 
(AFM) with load on the graphite mounted cantilever, 
very thin layers of graphite crystallites with lateral size 
~ 2 /im and thickness d raging from 10 - 100 nm are 
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FIG. 1: The inset shows an AFM image of a 12 nm thick 
mesoscopic graphite sample with four electrodes at the cor- 
ners for galvanomagnetic measurements. The left and right 
panel show the Hall resistance (R xy ) and magneto-resistance 
(R„) respectively as a function of magnetic field measured 
at T = 1.7 K in this device. Numbers near each curve indi- 
cate the applied gate voltages. In the right panel, curves are 
shifted for clarity and the dashed lines correspond to the zero 
lines of each curve. 



sheared off onto Si02/Si substrate. Multiple metal elec- 
trodes (Cr/Au) are then fabricated on the corners, using 
electron beam lithography (for AFM image of a typical 
device, see Fig. ^ inset). The degenerately doped silicon 
substrate serves as a gate electrode with thermally grown 
silicon oxide (500 nm) acting as the gate dielectric. 

Fig. d displays the Hall resistance (Rxy) and the 
MR (R X x) as a function of applied magnetic field, B, 
measured in a 12 nm thick graphite sample at tempera- 
ture T = 1.7 K. The excitation current is kept at 0.5 
/iA for both R xy and R xx measurements. The mag- 
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FIG. 2: Magneto-resistance measured in 12 nm (solid square) 
and 42 nm (open circle) thick sample at various gate volt- 
ages. The dashed and dotted lines are fits to a model de- 
scribed in text. The lower inset shows a schematic STB di- 
agram for electron (e e (fc)) and hole (e h (k)) in the presence 
of electrostatic potential induced by the gate. Upper inset 
represents the local magneto-resistivity across the sample, as- 
suming 0o = 8-Eo/M as an example. Symbols are defined in 
text. 



netic field is applied perpendicularly to the graphite basal 
plane. Both quantities exhibit oscillatory features on top 
of smooth backgrounds as B varies. Near V g R* V, 
the MR and Hall resistance exhibit similar behaviors to 
those observed in high quality bulk graphite The 
"V" shaped MR background is ascribed to the general 
nature of magnetotransport in materials with coexisting 
nearly compensated electron and hole carriers |14| , while 
the oscillations on top of the background are related to 
the SdH effect, the quantum oscillations due to Landau 
level formation Remarkably, as we vary the gate 

voltage, V g , the behavior of R xx and R xy changes dra- 
matically. The background in the MR is most prominent 
at V^ maa: « —15 V. As V g moves away from this value, the 
slope of the MR background becomes much smaller. The 
change of Hall measurement is even more drastic: R xy {B) 



changes its sign of overall slope as AV g 
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swings from negative to positive values, indicating that 
AV^ changes the dominant majority charge carriers from 
holes to electrons. This is a somewhat surprising result 
at first sight, since the thickness of the sample (12 nm) 
is still an order of magnitude larger than the screening 
length of graphite (A s « 0.4 nm (lfj), and thus only rela- 
tively small portion of the sample is affected by the gate 
electric field. We will discuss this point quantitatively 
below. 

The aforementioned EFE in mesoscopic graphite sam- 
ples is clearly presented by observing R xx as a function of 



gate voltage at a fixed magnetic field. Fig. [3 shows R xx 
as a function of V g at a large magnetic field (B m = 8.5 T) 
for two samples (d = 12 and 42 nm) As expected 

from Fig. Q Rxx has a peak near a gate voltage where 
AVg ~ 0, falling slowly as AV g \ becomes large. We found 
that this gate dependence strongly depends on d. For the 
12 nm sample, R xx is suppressed to ~ 10% of its peak 
value, while it is still ~ 60% for the 42 nm sample at 
AVg = 80 V. Such a sensitive dependence of R xx {V g )\B m 
on d is indicative of the reduced EFE by screening of 
induced charge near the sample surface. 

In order to elucidate the dependence of R xx on V„ , 
we employ the simple two band (STB) model |18j . 
which has been successful in understanding the MR in 
graphite 0, Q| . The STB model assumes that the bot- 
tom of the electron band and the top of the hole band 
overlap with a small band overlap 2Eq near the Fermi en- 
ergy Ep. The resistivity of a sample, p, in the presence 
of a magnetic field can be expressed by [lij : 



Ap _ 4[i 2 B 2 n e n h / {n e + n h ) 2 
Po 1 + \pB(n e - n/,)/ (n e + nh)} 2 



(1) 



where po = p(B = 0), Ap = p(B) — po, /i is the average 
carriers mobility, and n e and n/, are the carrier concen- 
trations of electrons and holes, respectively. Generally, 
Ap varies the most as a function of B when electrons 
and holes are nearly compensated (i.e., n e « n/,). From 
Fig. |21 we infer that this condition is met at V g ~ V™' 1 * 
where the growth of the MR background as a function 
of B is strongest in our samples (see the curves for 
V g = -10 V and Vg = -20 V in Fig. Q} 0. As AV g 
increases from zero, the induced charge in the sample 
screens the gate electric field and the electrostatic poten- 
tial in the sample is given by <f>(z) = 4>oe~ z f Xs , where z is 
measured from the interface between the sample and the 
substrate. The constant 4>o can be determined from the 
electrostatic gate coupling to the sample. By integrating 
over the induced charge in the sample, we obtain (f>o = 
aAVg with the constant or 1 = 1 + e (l - e~ rf/As )/ X s C g , 
where C g is the gate capacitance per unit area of the 
sample and £q is the vacuum permittivity 01 • 

We incorporate this local electrostatic potential to 
the STB model by considering a gradient in n e and 
71/,. Suppose AVg > 0, the local electrostatic potential 
will pull down the electron and hole bands by |e|0(z) 
(Fig- ID lower inset). For a sufficiently large gate volt- 
age, such that |e|0o > Eq, the sample can be devided 
into three regions by introducing a hole depletion depth, 
t dp = X s \og(\e\cj) Q /E ): (I) < z < t dp , where n h ~ 0; 
(II) tdp < z < tdp + X s , where < n/, < n e : and 



(III) z > tdp + A s , where 71/, 



In region (I), only 



electrons participate in the transport, and p(z) increases 
as z approaches zero, owing to the electric field induced 
accumulation of n e near the surface. In region (II), the 
MR is described by Eq. ^ so a steep increase of p{z) 
is expected as n e — nh becomes small. In region (III), 
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FIG. 3: (a) The SdH Oscillations observed in Fig-d after sub- 
traction of smooth backgrounds. Solid (open) symbols corre- 
spond to peak (valley) of the oscillations found after passing 
the curve through a low pass filter (dotted line). Curves are 
displaced for clarity, (b) Landau plots (see text) obtained 
from (a). Negative indices are assigned to the hole branch for 
clarity. Lines are linear fits to each set of points at different 
V g . Inset: the frequency of the SdH oscillations obtained from 
the slopes of the line fits in (b) as a function of gate voltage. 



the gate electric field is completely screened, so n e s=s 
and p(z) w p max = R xx (AV g = 0)d. The exact op- 
posite argument works for a sufficiently large negative 
gate voltage, where electrons are depleted. Note that 
for small |AVg|, where |e0o| < E Q , region (I) disap- 
pears (i.e., td p — 0). From above discussions, we now 
build a quantitative model to describe p{z). Accord- 
ing to the STB model n e (e), n^(e) oc e 3 / 2 , where e is 
measured from the bottom of the respective band edge, 
p(z)/pmax is obtained from Eq. ^ (Fig. [21 upper inset). 
Then the resistance of the sample can be evaluated from 
= p~ 1 (z)dz for fixed B and V g . Following in 
this way, a reasonable fit is obtained for both 12 nm 
(dashed line) and 42 nm (dotted line) samples as shown 
in Fig. |21 In this fit, we use Eq — 15 meV, a value quoted 
in [l|, and obtain C g = 26 aF//iin 2 (12 nm sample) and 
24 aF/^m 2 (42 nm sample) as a result of the fit. These 
capacitance values are in reasonable agreement with our 
previous estimations in a different analysis on the same 
samples 0. It is noteworthy that for a large V g such 
that \e\4>o » E , p(z) -C p max in region (I), and thus 
a significant portion of the total current flows in this 
region. Furthermore, as AV^ increases, td P grows only 
logarithmically. Even at AV g w 100 V, the largest gate 
voltage applied, td P ~ 1 nm, which corresponds to only 
~ 3 bottom layers. Therefore, only a few of the bottom 
layers of the sample are responsible for the observed EFE 
modulation of the galvanomagnetic transport quantities. 



We now turn our attention to the quantum oscillations 
observed in our mesoscopic graphite samples. The strong 
EFE modulation of the carrier density in the bottom lay- 
ers allows us to probe the quantum oscillations in these 
layers with a continuously tunable carrier concentration. 
Fig. |3f a) redisplays the separated SdH oscillations as a 
function of B^ 1 , obtained from the MR data shown in 
Fig. U after subtracting out the smooth background. The 
SdH oscillations indicate the oscillatory density of states 
at Ef as a quantized Landau level passes through Ef- 
The frequency of SdH oscillations, f s , is related to the ex- 
tremal area of the electron and hole pockets of the Fermi 
surface by f s — ficA e ^ h /2ir\e\, where A^. and A\ are the 
areas of extremal electron and hole pockets, and H and 
c are Plank constant and speed of light respectively 0] . 
Since n e and rih are modulated by V g , the observed vari- 
ation of f s can be explained by the change of A e ^ h . 

In order to demonstrate the change of f s quantita- 
tively, we first locate the major peaks (solid symbols) 
and valleys (open symbols) in the SdH oscillations after 
low pass filtering of the data jl^. The value of B~ x for 
a peak (valley), -B" 1 , is indexed by v, an integer (a half 
integer) number that corresponds to the Landau level re- 
sponsible for the particular oscillation. Fig. 02b) shows 
that each set of points (B^,v) at a given V g are on a 
straight line that intercepts the origin, implying that the 
period of SdH oscillations is regular. From the slope of 
these lines we obtain f s at different V g (Fig.|3fb) inset). 
The obtained / s 's are increasing with |AVg|. Therefore, 
we believe that the obtained f s corresponds to At for 
AV g > and to A\ for AV g < 0. This conclusion allows 
us to compare the experimentally observed f s with the 
expected value from the STB model. Assuming the Fermi 
surface of graphite is described by the overlap of elec- 
tron and hole bands in STB model, A\ ex (a AVg + Eq) 
and A\ oc (—aAV g + E ). This relationship leads to 
f s (\AV g \)/f° = l+a\AV g \/E 0i where /° = f s (AV g = 0). 
From the values of a and Eq, determined separately 
above, we estimate f{AV g = 50V)/ fo w 4.8, which is in 
reasonable agreement with the experimentally observed 
ratio 4.3. 

Finally, we discuss the temperature dependence of the 
SdH oscillations. Fig. 0] shows the oscillatory MR at two 
extreme gate voltages, V g — +40 V and V g = —60 V, 
at various temperatures. At these extreme gate voltages, 
the transport in the sample is dominated by only one 
type of carriers in a few bottom layers. Thus, the SdH 
oscillations in the upper (lower) panel of the figure cor- 
respond to electron (hole) Landau levels in the sample. 
In both cases, the observed SdH oscillation amplitude 
is gradually damped away as the temperature increases. 
The temperature dependent SdH oscillation amplitude 
has been used to extract the effective mass of charge car- 
riers [2]]]. At a fixed magnetic field, the temperature 
damping factor of the SdH oscillation amplitude is given 
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FIG. 4: Normalized magneto-resistance of the sample in FigQ 
at V g = 40 V (upper panel) and V g = — 60 V (lower panel). 
For the upper panel data are taken at temperatures 1.7, 4, 
10, 13, 15, 20, 25, and 30 K. For the lower panel, data are 
taken at temperatures 1.7, 4, 10, 13, 15, and 20 K. Insets: 
SdH oscillation amplitudes divided by temperature, T, at a 
fixed magnetic field at above temperatures. The solid lines 
are fits to a model (see text). 

by: 

= 2K 2 k B Tm*/ehB 
T sinh(27r 2 k B Tm*/e%B) [ ' 

where to* is the effective mass of the carriers. We find 
that Eq.|21fits the observed amplitude damping very well 
(Fig. 01 insets). As a result from the fittings, the effective 
electron mass to* = (0.052 ± 0.002)TO e and hole mass 
m h = (0-038 ± 0.002)to c are obtained, where m e is the 
bare electron mass. These values agree well with 0.057TO e 
and 0.039to £ , reported effective mass of electrons and 
holes in high quality bulk graphite crystal [lj. 

In summary, we report galvanomagnetic transport 
in mcsoscopic graphite samples consisting of tens of 
graphene layers. Strong modulation of the Hall resistance 
as well as the magneto-resistance has been observed as 
the applied gate voltage changes. The Landau level for- 
mation of electron and hole carriers is also tuned by the 
gate. The unique experimental method discussed here 
can be applied to other layered materials to investigate 
novel transport phenomena in unconventional two dimen- 
sional systems. 
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